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I.  ELECTRICAL  AND  OPTICAL  PROPERTIES 

INFRARED  MEASUREMENT  OF  HOT  CARRIER  TEMPERATURE  IN  p-TYPE  GE 

W.  Pinson  and  R,  Bray 

Since  the  last  report  a  review  and  analysis  of  the  methods  and  results  of  the 
infrared-hot  carrier  experiment  have  been  completed,  which  indicated  the  desirability 
of  making  many  modifications  in  the  apparatus  and  improvements  in  the  procedure  for 
taking  the  data  and  analysing  it.  Much  of  the  new  equipment  has  now  been  designed, 
built  and  tested.  The  highlights  of  the  above  work  will  be  outlined  below. 

The  change  in  absorption  coefficient,  Act,  with  pulsed  electric  field  strength 
will  be  directly  determined  by  measuring  the  ratio  of  the  pulse  change  in  transmis¬ 
sion,  A  I.J.,  due  to  a  pulsed  field,  to  the  transmission,  1^,  at  zero  field.  The 
expression  relating  A  a  to  this  ratio  is  a  relatively  simple  one  and  contains  all  the 
necessary  corrections  for  multiple  reflections  and  the  fact  the  light  path  through 
the  sample  is  greater  than  the  region  subjected  to  high  fields.  The  latter  geomet¬ 
rical  consideration  is  due  to  the  presence  of  large  area  end  contacts  to  eliminate 
injection  effects  at  high  fields.  The  new  procedures  also  eliminate  possible  errors 
due  to  long  time  variations  in  light  intensity  and  detector  sensitivity. 

However,  the  ratio  A  Ij/lj  cannot  be  measured  directly  because  of  1)  the  non¬ 
linearity  of  the  detector,  2)  the  long  excited  carrier  lifetime  t  of  the  InSb 
detector  (t  >  pulse  length),  3)  lattice  heating  effects  produced  by  the  pulse  in  the 
Ge  sample.  Improved  techniques  and  analyses  for  all  these  effects  have  been  devel¬ 
oped. 

Our  determination  of  the  absorption  coefficient  a  as  a  function  of  wavelength, 
particularly  at  77^K,  and  the  changes  in  (X  with  lattice  temperature  had  not  been  too 
successful  for  Ge  with  around  3  x  10^^  holes/cc  in  the  2  to  5  micron  region.  The 
main  source  of  the  difficulty  in  measurement  lies  in  the  small  values  of  a  in  this 
range  (0.01  cm  ^  to  1 .0  cm  ^),  necessitating  the  use  of  long  samples  with  the  atten¬ 
dant  problems  of  uniform  cooling  and  optical  defocussing.  We  have,  therefore,  built 
an  exchange  gas  cell  to  cool  the  samples,  and  switched  to  a  technique  of  making  com¬ 
parative  transmission  measurements  of  our  p-type  sample  with  a  pure  n-type  Ge  sample 

of  exactly  the  same  length,  resulting  in  an  improvement  of  our  measurements  of  a.  As 

I  2 

the  several  published  ’  values  of  the  magnitude  of  cc  disagree  by  30^  or  more,  it  is 
necessary  for  us  to  remeasure  absorption  cross  sections  more  carefully. 
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1.  H.  B.  Briggs  and  R.  C.  Fletcher,  Phys.  Rev.  9J.»  (1953). 

2.  R.  Newman  and  W.  W.  Tyler,  Phys.  Rev.  105.  885  (1957). 


GALVANOMAGNETIC  PROPERTIES  OF  n-GaP 
W.  M.  Becker  and  H.  Y.  Fan 

Optical  measurements  by  Edwards  and  Drickamer^  Indicate  that  the  lowest  conduc¬ 
tion  band  in  n-GaP  is  silicon-like,  with  the  band  minima  lying  along  [100]  directions. 
However,  no  extensive  electrical  measurements  are  available  relating  to  this  result. 
We  have  begun  an  investigation  into  the  gal vanomagnet ic  properties  of  this  material 
using  sulphur-doped  single  crystals  supplied  by  Bell  Telephone  Laboratories. 

initial  work  showed  that  the  usual  solder  or  pressure  contacts  resulted  In 
unsatisfactory  non-ohmic  high  resistance  contacts.  Successful  contact  to  the  samples 
was  finally  made  as  follows:  The  appropriate  regions  of  the  sample  surface  were 
painted  with  Hanovla  #38  silver  paste  to  which  a  small  amount  of  tellurium  had  been 
added.  The  sample  was  then  heated  in  a  hydrogen  atmosphere  at  600^C  for  one  hour, 
after  which  the  silver  paste  was  removed.  Low  resistance,  ohmic  contacts  were 
obtained  by  soldering  indium  directly  to  the  diffused  areas. 

Hall  effect,  conductivity,  and  magnetoresistance  were  measured  on  a  sample 

3  2 

having  a  room  temperature  Hall  coefficient  R  =  -19  cm  /coul .  and  Pj,  =  131  cm  /volt- 

A  o  n 

sec.  At  77^K,  the  Hall  coefficient  rises  to  R  =  -2.2  x  10  cm  /coul  and  =  333 

2  ” 
cm  /volt“Sec.  Transverse  magnetoresistance  was  observed  neither  at  room  temperature 

for  fields  up  to  22,000  gauss,  nor  at  77^K  for  fields  up  to  13,500  gauss.  Based  on 
the  errors  in  the  measurements,  the  room  temperature  value  was  estimated  to  be  of  the 
order  of  Ap/p  ^  O.ljK  at  the  highest  field;  the  estimated  value  at  77^K  was  Ap/p  ^  \ 
0,55^.  Measurements  at  higher  fields  should  yield  more  reliable  data  for  the  deter¬ 
mination  of  the  magnetoresistance  coefficients. 


1.  A,  L,  Edwards,  T.  E.  Slykhouse,  and  H.  G.  Drickamer,  J.  Phys.  Chem.  Solids  n_,  140 
(1959). 
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INFRARED  ABSORPTION  IN  GaSb 
E.  J.  Johnson  and  H.  Y.  Fan 

Three  exclton  like  lines  observed  on  the  fundamental  absorption  edge  of  GaSb  at 

4,5°K  have  been  described  In  previous  progress  reports  These  lines  have  now 

been  observed  at  1.7°K  where  they  become  sharper  and  more  intense.  The  lines  at 

1.7°K  are  shown  in  Fig.  1.  The  a  line  is  probably  associated  with  the  formation  of  a 

free  exciton.  A  different  interpretation  is  necessary  for  the  p  and  y  lines.  The 

possibility  under  consideration  is  that  the  latter  lines  might  be  due  to  the  forma- 

2 

tlon  of  excitons  near  imperfections  .  The  obvious  test  for  this  interpretation  would 
be  to  look  for  variation  in  the  strengths  of  the  lines  in  different  samples.  In  an 
attempt  to  verify  this  hypothesis,  measurement  was  made  on  a  sample  having  a  somewhat 
smaller  Hall  coefficient.  The  result  showed  that  the  y  peak  broadened  into  a  step 

without  appreciably  increasing  the  strength.  It  is  possible  that  the  excess  impurity 

in  the  sample  is  not  of  the  type  that  gives  rise  to  the  y  peak.  Preliminary  measure¬ 
ments  of  the  Zeeman  effect  have  been  made  on  the  p  and  y  lines.  The  results  obtained 

are  shown  in  Fig.  2  and  can  be  summarized  as  follows: 

1.  Each  line  splits  into  two  components. 

2.  Each  component  occurs  at  higher  photon  energy  than  the  line  at  zero  magnetic 
field. 

3.  The  spectra  with  polarization  E)|H  is  identical  to  the  spectra  with  EXH. 

4.  The  observed  splitting  is  too  large  to  be  accounted  for  by  spin. 

These  results  need  further  study,  but  they  seem  to  indicate  effects  of  degen¬ 
erate  valence  bands  at  k  =  0. 


1.  E..  J.  Johnson  and  H.  Y.  Fan,  Semiconductor  Research,  Third  Quarterly  Report, 
April  1,  1961  to  June  30,  1961. 

2.  E.  J.  Johnson  and  H.  Y.  Fan,  Semiconductor  Research,  Final  Report,  October  1, 
1960  to  September  30,  1961. 
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RECOMBINATION  RADIATION  IN  GALLIUM  ANT  I  MON  IDE 
I.  FlllnskI  and  H,  Y,  Fan 

Investigation  of  recombination  radiation  in  GaSb  has  been  continued.  In  the 
previous  measurements  at  liquid  helium  temperature,  there  was  some  uncertainty 
regarding  sample  temperatures.  The  measurements  were,  therefore,  repeated  both  at 
liquid  nitrogen  and  liquid  helium  temperatures  by  actually  immersing  the  sample  in 
the  coolant.  The  present  results  show  that  the  position  of  peaks  observed  at  liquid 
helium  temperature  for  pure  p-type  samples  are  0.752  ev  +  0.002  ev,  0.779  +  0,02  ev, 
and  0.796  +  0.02  ev.  The  half  width  of  observed  peaks  is  ^  0,006  -  0.008  ev  and  is 
larger  than  the  resolution  used  viz  ^  0.002  ev.  The  relative  intensity  of  0.779  ev 
to  0.796  ev  peaks  changed  from  sample  to  sample  even  for  samples  taken  from  different 
parts  of  the  same  ingot. 

The  following  model  has  been  considered  to  explain  the  observed  emission  and 
absorption  spectra.  The  0.779  ev  and  0.752  ev  peaks  are  connected  with  electron 
recombination  with  acceptor  centers.  The  0.796  ev  peak  is  interpreted  as  exciton 
recombination  at  some  impurity.  This  should  be  compared  with  the  free  exciton 
recombination  expected  to  occur  in  purest  samples  with  energy  ~  0.811  ev  on  the  basis 
of  absorption  measurements  of  E.  J.  Johnson  and  H.  Y.  Fan  reported  in  this  progress 
report. 

MAGNETORESISTANCE  IN  N-TYPE  GERMANIUM  IN  THE  IMPURITY  CONDUCTION  RANGE 

G.  Sadasiv 

The  transverse  magnetoresistance  of  several  n-type  Germanium  samples  has  been 
measured  in  the  temperature  range  between  4,2^K  and  1 .2^K,  with  magnetic  fields  up  to 
12,000  oersted.  The  orientation  of  the  samples  and  the  characteristics  at  room 
temperature  and  at  77.4^K  are  given  in  Table  I. 
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TABLE  I 


297° 

K 

77.4° 

K 

Sample 

Impurity 

<i> 

<  H  > 

3  ''  ,-l 

cm.  coul . 

P 

ohm-cm. 

3  -1 

cm.  coul . 

P 

ohm-cm 

1 

Sb 

001 

1 10 

10.5 

. 00848 

13.2 

.0109 

2 

tl 

001 

1 10 

44.3 

.0254 

70.8 

.0247 

3 

II 

on 

100 

57.4 

.0296 

92.1 

.0251 

4a 

tl 

001 

110 

76.9 

.0382 

119.7 

.0293 

5 

tl 

001 

1 10 

90.9 

.0410 

137.5 

.0292 

6 

tl 

001 

1 10 

106.8 

,0480 

168.2 

.0305 

7 

.11 

001 

1 10 

,0687 

.0367 

8 

As 

on 

oTi 

18.9 

.0144 

29.3 

.0204 

9 

tl 

on 

oTi 

25.8 

.0167 

41.4 

.0219 

For  sariiples  5  through  9,  measurements  were  also  made  with  the  magnetic  field  along 
the  <  100  >  direction.  There  was  no  significant  difference  in  the  magnetoresistance. 
For  samples  3,  4a,  5,  and  9  the  resistivity  in  a  magnetic  field  H,  at 
,  temperatures  between  4.2°K  and  1 .2  n,  could  be  described  by  the  relation 

Ph  =  pQ  ®xp  (e  +  a  H^)/kT  (1) 

For  sample  2,  the  same  type  of  relation  held  for  fields  greater  than  8000  oersted. 
Sample  6  showed  the  62  anomaly  in  the  resistivity^  and  in  the  62  region  the  magneto- 
resistance  was  described  by  (1).  The  values  of  Pq,  c,  and  a  are  given  in  Table  M. 

TABLE  II 


Sample 

PO 

ohm-cm. 

e 

ev 

“  -2 
ev  oersted 

2 

.143 

.030  X  10"^ 

-13 

.35  X  10 

3 

.297 

.113 

1 .27 

4a 

.532 

.  .284 

2.35 

5 

.748 

.691 

3.53 

6 

.698 

1 .25 

4.05 

9 

.223 

.635 

1.81 

6 
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In  Fig.  3,  the  values  of  Ae  (=  a  H  )  are  plotted  against  H  for  these  samples. 

Sample  1  had  a  negative  magnetoresistance,  reported  in  a  previous  progress 
2 

report  .  For  sample  7  A  p/p  at  a  constant  field  decreased  as  the  temperature  was 
lowered,  a  behaviour  similar  to  that  observed  in  pure  samples  by  Fritzsche*. 

1.  H.  Fritzsche,  Phys.  Rev.  99,  406  (1955). 

2.  G,  Sadasiv,  Semiconductor  Research,  Third  Quarterly  Report,  April  I  to  June  30, 
1961 . 


CYCLOTRON  RES0NANCE-5mm  APPARATUS 
T.  Tohver  and  G.  Newell 

Some  difficulties  in  the  design  and  construction  of  5  mm  sample  cavities  have 
been  overcome.  Initially,  to  facilitate  drilling  coupling  holes  and  machining 
sample  illumination  slits,  the  cavity  had  been  split  into  two  parts.  However,  at 
5  mm  wavelengths,  machining  imperfections  made  this  impractical;  the  split  interfered 
with  the  wall  currents  and  power  was  lost  through  the  Illumination  window,  both 
effects  combining  to  lower  the  Q  excessively.  An  improved  design  eliminates  the 
split  and  the  illumination  window  consists  of  a  series  of  holes  beyond  cutoff  (Fig.  4). 
The  new  design  has  been  tested;  it  gave  a  cavity  Q  of  ~  2000  which  is  quite  satis¬ 
factory.  The  variable  coupling  scheme  shown  on  the  accompanying  sketch  (Fig.  4)  also 
works  quite  well  -  varying  smoothly  from  zero  coupling  to  overcoupling  between  the 
waveguide  and  the  cavity.  These  results  on  the  sample  cavity  and  the  latest  perfor¬ 
mance  data  on  5  mm  video  detectors  have  permitted  a  rough  estimate  of  the  sensitivity 
of  the  equipment.  These  calculations  are  included.  Also  included  is  a  schematic 
diagram  of  the  microwave  circuit  with  an  indication  of  the  sections  that  have  been 
completed  (Fig.  5). 

Sensitivity  Estimate 

Optical  excitation  of  the  carriers  is  to  be  used.  The  sensitivity  of  the  appara¬ 
tus  will  be  defined  as  the  minimum  number  of  detectable  carriers.  Straightforward 
circuit  theory  gives  for  the  detected  power  under  optimum  cavity  matching  conditions 

P .  -  Po  X  0.193 

a  up 


Pj  -  detected  power 

Po  -  power  incident  on  sample  cavity 

AQ  -  change  in  cavity.il  due  to  excited  carriers 

Qo  unloaded  cavity  Q,. 

Assuming  the  power  lost  to  the  excited  carriers  is  small  compared  to  the  wall  losses 


Qo 


Sample  „ 

as  "  °  ^  '' 

2)0  Cavity  « 

;  .  d  V 


a  -  conductivity  due  to  excited  carriers  only. 

For  a  needle-shaped  sample  with  a  diameter  small  compared  to  the  wavelength  and 
placed  parallel  to  the  maximum  E-fleld  of  a  cavity 

Sample  „  , 

r  a  d  V  Kdc 

max  Samp  I  e 


and 


Cavity 


d  V  rw  1/4  * 

o 


E  V 
max  Cav i ty 


^  »  222_  ^Sample 

Qo  to  e  V*  , 

^  o  Cavity 


Now  a  V-  ,  *  q  (p  Ap  +  u  An)  V_ 

Samp I e  ^  P  e  S 

For  an  order  of  magnitude  calculation  assume 


/  2 
,  «4  cm 

u  u  “10  — r-—— 

e  p  volt-sec 

2 

.  I  - £0 - 

vol t-sec 

Also  cal  1 


Vg  X  (Ap  +  An) 


N 

Total  number  of  carriers  excited. 
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m  “  3.3  X  10'*  rad/sec 
“  8.9  X  10  farad/m 

q  “  1.6x10^^  coul 

'^Cavity  ~  (mm)^  *=  10"^ 

Po  ~  1  mill iwatt 

Then  ~  2  x  10"’^  N  watts. 

At  this  point,  the  detector  characteristics  enter. 

2 

A  crude  equivalent  circuit  for  the  square  law  detector  is  given  below 

3 

Typically  P  «  0.13  amps/watt 

R  -  15  K 

o 

Then  the  signal  voltage  is  given  by: 

Vj  -  1$  Rq  “  P  ’’d  ^  10”*^  N  x  0.13  X  15  X  10^. 

Thus  Vg  »- 4  X  10"^®  N  volts 

Assume  the  1/f  noise  voltage  across  R^  at  330  c.p.s.  is  100  times  the  thermal  noise 
2  ° 

in  R  . 
o 

Hence 


“  100  X  (K  T  R^ 
noise  o 

For  a  30  cps  bandwidth  (AF) 

V  ~  4  X  lO”®  volts 
n 

For  the  minimum  detectable  signal 

y.min  ^  y 

signal  noise 

Vj*"  »  4  X  10“^  V  -  4  X  10"’ ° 

^min 


4 

10  carriers 
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1.  G.  Ffeher,  B.S.T.J.  36,  449  (1957). 

2.  E.  L.  GInzton,  Microwave  Measurements  (McGraw-Hill  1957),  p.  117-119, 

3.  Sylvania  04074  Millimeter  Wave  Detector  (Private  Communication) (B.  R.  Lincoln, 
Produce  Sales  Specialist). 
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11.  IRRADIATION  EFFECTS 

ELECTRON  PARAMAGNETIC  RESONANCE  -  SILICON 
W.  Jung  and  G.  Newell 

A  sample  of  neutron  irradiated  crucible  grown  silicon  (sample  #2^  USSC  No.  4, 

n  1  Q 

100  Cl-cm;  ORNL  Irradiation  at  T  ~  50  C,  0  =*  1 .4  x  10  nvt)  has  been  carried  through 
an  isochronal  anneal  in  the  following  steps:  90°  to  170°C  in  10°  steps  for  30 
minutes;  170°  to  185°C  in  5°  steps  for  30  minutes;  and  200°  to  500°C  in  50°  steps  for 
60  minutes .each.  These  heat  treatments  were  carried  out  in  a  quartz  furnace  in  a 
helium  atmosphere. 

After  each  step  of  annealing,  a  complete  measurement  of  the  angular  variation  of 
the  EPR  spectrum  at  room  temperature  was  made,  using  intervals  of  2.5°  for  critical 
runs  and  10°  otherwise.  The  microwave  Q  remained  high  through  steps  0  to  17  (400°C) , 
indicating  that  the  Fermi  level  remained  in  the  middle  of  the  gap,  as  expected  from 
previous  work.  The  Q,  decreased  somewhat  after  step  18  (450°C) ,  and  quite  appreciably 
after  step  19  (500°C),  at  which  point  only  unanalysable  traces  of  resonance  remained. 
The  centers  formed  at  or  below  400°C  are  thus  characteristic  of  intrinsic  material, 
but  the  Fermi  level  after  the  450°  anneal  is  uncertain.  A  measurement  of  its  final 
Fermi  level  is  under  way. 

Growth  and  Decay  Curve: 

Some  12  to  14  centers  are  observed  to  grow  and  decay  during  annealing.  Fig.  6 
plots  the  concentration  of  those  centers  for  which  g-tensors  have  been  worked  out  and 
which  can  therefore  be  positively  identified.  As  discussed  below,  the  curves  for  a 
single  center  are  internally  consistent  to  +  \5^,  while  the  relative  and  absolute 
concentrations  are  uncertain  by  factors  of  2  or  3  either  way. 

These  curves  clarify  the  relations  among  the  previously  identified  centers  (N, 

I -III)  and  also  the  relation  between  Oak  Ridge  (T  50°C)  and  Argonne  (T  100°C) 
irradiated  samples.  Center  I  was  originally  observed  in  Argonne  irradiated  material 
containing  arsenic  (Sample  #1 S ;  it  is  now  found  in  pure  crucible  grown  material 
irradiated  at  the  lower  temperature  and  annealed  at  300  to  350°C.  In  fact,  having 
determined  the  g-tensors  for  centers  I',  V,  and  VI,  we  have  now  identified  all  the 
important  lines  in  the  Argonne  sample  and  find  that  their  relative  intensities  cor¬ 
respond  roughly  to  the  300°C  step  of  annealing. 

A  striking  feature  of  the  annealing  curves  shown  In  Fig.  6  is  the  occurence  of 
pairs  of  centers  which  grow  and  decay  together,  and  whose  g-tensors  show  marked  simi- 
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lari  ties,  as  discussed  later.  The  pairing  of  centers  V  and  VI,  for  example,  can 
hardly  be  coincidental,  but  no  detailed  explanation  is  yet  available. 

In  these  complicated  spectra,  lines  overlap  more  often  than  not.  Using  the 
angular  variation  plots  for  each  run,  all  the  lines  which  appear  (accidentally)  to 
be  in  the  clear  at  each  angle  are  picked  out  for  measurement,  giving  15  and  40 
intensity  measurements  for  each  center  and  temperature.  These  data  are  used  to  con¬ 
firm  the  g-tensor  assignments  and  to  estimate  defect  concentrations.  The  accuracy  of 
these  estimates  is  discussed  below. 

(a)  Confirmation  of  g-tensors:  Each  <  110  >  type  center  is  represented  by  seven 
branches  in  the  angular  variation  plot.  To  verify  that  the  branches  assigned  to  a 
given  center  do  in  fact  grow  and  decay  together,  separate  annealing  curves  were  con¬ 
structed  for  each  branch.  For  a  given  run  the  branches  belonging  to  a  single  center 
showed  a  maximum  spread  of  2:1  in  intensity,  with  an  r.m.s.  deviation  of  ]6^  from 
the  mean.  As  several  of  the  centers  were  followed  over  a  20:1  variation  in  concentra¬ 
tion,  the  agreement  is  convincing.  Taken  with  the  very  precise  internal  consistency 
criteria,  this  seems  sufficient  to  confirm  the  g-tensor  assignments.  However,  the 
possibility  that  pairs  of  g-tensors  belong  to  a  single  center  is  mentioned  below. 

(b)  Relative  concentrations:  The  standard  deviation  of  the  mean  intensity 
calculated  for  each  center  is  +  5  to  10^6,  and  is  less  than  the  systematic  variation 
of  +  15^  seen  in  the  curves  for  centers  II  and  III  below  130^C.  This  variation 
appears  to  result  from  small  relative  changes  in  the  filling  factors  of  standard  and 
sample  when  the  latter  Is  replaced  after  each  run.  The  annealing  curve  for  a  given 
center  is  thus  self  consistent  to  +  15  or  20^^. 

Assuming  similar  line  shapes  (gaussian),  the  relative  concentrations  of  two 
centers  are  given  by: 

Hi  ^  AnA' 

"2  ■  *2’’  {^"ij  • 

where  A  Is  the  amplitude  and  A  H  the  line  width.  The  20^  uncertainty  in  A  H  intro¬ 
duces  an  uncertainty  of  1.4:1  in  the  concentrations.  Further,  the  amplitude  of  the 
hyperfine  satellites  must  be  added  to  that  of  the  principal  line.  This  correction 
amounts  to  +  50^6  for  Center  II,  but  has  not  been  applied  as  the  h.f.s.  of  the  other 
centers  is  as  yet  unknown.  Relative  concentrations  are  thus  reliable  only  within  a 
ratio  of  2:1  unless  the  hyperfine  structure  can  be  assumed  similar. 

(c)  Absolute  concentrations:  The  absolute  concentration  scale  is  derived  by 
comparison  with  the  internal  ruby  standard  using  the  relation 
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Ah2Y  ^ 
^  ‘  V^y  V2 

-  !  .24  X  10’^  (A  HJ^  -r 


The  subscripts  1  and  2  refer  to  ruby  and  the  sample,  respectively,  while  the  various 
quantities  are  defined  as  follows: 


!1l 

% 


<  "rf  >1 

<  >2 


0.216 


1 

(— )  ■»  transition  probability  ratio  ~  4 

“2 

A  H2  “  line  width  of  ruby  *  15.9  gauss 

1 6 

Nj  *  number  of  Cr  spins  In  the  Ruby  standard  ^2.1  x  10  . 


•2 

=  sample  volume  •  4,12  x  10  cc. 


“  *  amplitude  ratio 

1 


1 1 1 

Nj  is  estimated  from  the  nominal  O.ljJCr  concentration  in  the  ruby  and  may  be  In 
error  by  a  factor  of  2.  Until  a  calibration  of  the  standard  is  completed,  the  abso¬ 
lute  concentrations  are  uncertain  to  a  factor  of  3  or  4:1. 
g-tensors 

The  g-tensors  and  other  properties  of  the  centers  discussed  above  are  collected 

2 

In  Table  1.  The  tabulated  values  replace  those  given  In  the  last  report  ;  in  particu¬ 
lar,  improved  values  for  Center  I  have  been  obtained  from  the  new  data,  while  the  N 
center  g-tensor  has  been  converted  to  absolute  values  consistent  with  the  rest  of  the 
table*  Center  IV  has  not  been  identified  with  certainty  in  this  sample,  but  is 
included  for  completeness.  With  the  exception  of  IV,  VII,  and  VIM,  all  the  centers 
listed  have  been  observed  in  at  least  two  independent  samples.  Centers  II,  III,  and 
N  appear  In  floating  zone  material;  the  oxygen  dependence  of  the  others  is  unknown 
at  present. 

To  obtain  these  g-tensors,  a  cos  *12(0  -  ]“  curve  is  fitted  to  each  branch  of 

the  angular  variation  plot  to  determine  adjusted  values  for  the  nine  different  g- 
values  belonging  to  the  three  principal  directions.  Four  of  these  values  are  used  to 
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compute  the  g-tensor;  the  remainder  are  subjected  to  consistency  checks.  The  result¬ 
ing  g-tensors,  in  general,  fit  the  data  to  +  0.5  gauss,  or  within  half  a  line  width. 

A  simplified  least  squares  treatment  is  being  performed  to  improve  the  g-tensors;  the 
adjustments  are  not  expected  to  exceed  the  quoted  uncertainties  of  .0002  in  g  and 
0.2*^  in  0.  The  g-value  scale  is  derived  from  simultaneous  measurements  (with  the  same 

counter)  of  the  microwave  frequency  and  the  proton  precession  frequency,  and  is  free 

3 

of  systematic  error  other  than  that  introduced  by  our  use  of  the  published  ratio 

\i 

»  657.4758, 

which  is  valid  for  protons  in  mineral  oil,  whereas  our  proton  probe  uses  glycerine. 
Discussion: 

Fig.  1  shows  a  striking  tendency  of  the  damage  centers  to  occur  in  pairs  (II 
and  III,  V  and  VI,  I  and  I*,  VII  and  VIII)  which  grow  and  anneal  together.  Examina¬ 
tion  of  the  tabulated  g-tensors  shows  a  general  similarity  between  members  of  a  pair; 
this  similarity  is  rather  more  striking  when  the  angular  variation  plots  are  examined. 
As  this  report  is  written  we  speculate  that  either  the  pairs  represent  two  configura¬ 
tions  of  the  same  basic  center,  or  that  the  two  g-tensors  actually  represent  the 
m^  =  -1  =  0  and  m^  =  0  -*  m^  =*  +1  transitions  of  the  same  S  «  1  system.  We  have 

not  yet  had  time  to  test  these  possibilities. 

In  addition  to  the  centers  tabulated,  at  least  two  more  are  present  in  the 
annealing  stages  near  170^C,  and  one  or  two  more  in  the  450^C  stage.  They  are  heavily 
obscured  and  it  remains  to  be  seen  whether  further  analysis  will  disentangle  them. 

Several  questions  for  future  Investigation  are  raised  by  these  results.  Fore¬ 
most,  perhaps,  is  the  determination  of  the  hyperfine  structures  In  order  to  obtain 
information  bearing  on  the  microscopic  structures.  Also,  however,  there  is  the 
question  of  oxygen  dependence;  the  rapid  and  complex  phenomena  in  the  450^C  region 
suggest  a  connection  with  the  coagulation  of  oxygen  impurity  and  require  more  detailed 
investigation.  We  have,  however,  abandoned  the  investigation  of  reaction  kinetics 
suggested  in  the  last  report:  the  system  appears  far  too  complex* 
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TABLE  I 


Center 

~  T°C 

~  T°C 

(a) 

^1 

92 

93 

e 

AH 

Impurity 

forms 

anneals 

i.0002 

+.0002 

+.0002 

+0.2° 

+0.2  g. 

dependence 

1 1 

? 

180° 

2.0109 

2.0127 

1 .9974 

14.9° 

B 

i ndependent 

1 1 1 

? 

180° 

2.0010 

2.0072 

2.0128 

32° 

B 

II 

N 

150° 

500° 

2.0089 

2.0125 

2.0045 

0 

00 

1 .2 

II 

1 

250° 

400° 

2.0194 

2.0243 

1 .9785 

12.8° 

B 

i ndependent 

1' 

250° 

400° 

1 .9971 

1  .9942 

2.0305 

2.8° 

B 

or 

Oxygen 

V 

250° 

450° 

2.0125 

2.0139 

1 .9940 

4.7° 

1 .2 

11 

VI 

250° 

450° 

2.0015 

2.0041 

2.0148 

8.0° 

1 .2 

II 

VII 

300° 

500° 

2.0127 

2.0149 

1 .9937 

0° 

1.6 

It 

VI 1  1 

300° 

500° 

2.0017 

2.0031 

2.0167 

0° 

1.6 

II 

IV 

400°? 

500° 

2.0088 

2.0113 

2.0050 

22.0° 

1.1 

It 

(a)  02  is  along  <  ITO  >;  lies  in  the  plane  containing  <001  >  and  <  110  >, 
making  an  angle  6  with  <  OOi  >. 

(b)  AH  =  full  width  between  inflections,  in  gauss. 


1.  The  samples  used  in  this  work  are  described  in  a  previous  report:  Jung  and 
Newell,  Semiconductor  Research,  3rd  quarterly  report,  April  1  to  June  30,  1961, 
p.  13. 

2.  Jung  and  Newell,  Semiconductor  Research,  Final  Report,  Oct.  1,  1960  to  Sept.  30, 
1961,  p.  10. 

Cohen,  Du  Mont,  et  al..  Rev.  Mod.  Phys.  27,  363  (1955). 


3. 
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ELECTION  PARAMAGNETIC  RESONANCE  -  GERMANIUM 
D.  Trueblood  and  G.  Newell 

I .  Apparatus 

Construction  of  the  x-band  superheterodyne  spectrometer  is  continuing.  The 
microwave  circuit  and  the  klystron  control  circuits  are. complete  and  operate  satis¬ 
factorily.  Amplifiers  and  demodulators  for  the  signal  channel  have  been  completed 
and  are  undergoing  test, 
y .  Filling  Factor  Calculations; 

As  mentioned  in  the  last  report,  measurements  will  be  made  on  cylindrical 
germanium  samples  located  on  the  axis  of  a  silicon  filled  TE^jj  mode  cylindrical 
cavity.  The  configuration  Is  sketched  In  Fig.  7,  which  also  gives  calculated  curves 
of  filling  factor  and  cavity  Q  as  a  function  of  sample  diameter.  These  curves  serve 
to  determine  the  optimum  diameter  for  a  sample  of  given  conductivity.  They  supple¬ 
ment  the  results  previously  given  for  rectangular  cavities'.  The  derivation  of  the 
expressions  valid  for  small  sample  diameters  is  sketched  below. 

The  rf  field  configuration  In  the  cavity  for  r  ^a  is  given  by: 
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where  a  is  the  radius  of  the  Ge  sample,  Xqj  “  3.832  for  a  TEqj  mode,  is  diameter 

of  an  air  filled  cavity,  k  is  the  dielectric  constant,  is  an  amplitude  factor, 

X  Is  free  space  wavelength,  X  is  the  guide  wavelength,  and  D  the  loaded  cavity 
o  g 

diameter. 

A  and  D  were  calculated  as  functions  of  a  by  matching  .  and  H  across  the 
o  in  tan. 

boundary  r  *  a.  X  also  depends  on  a,  but  for  small  values  of  a  both  D  and  X. 

9  9 

approximately  equal  the  values  for  a  Si  filled  cavity. 

The  filling  factor  integral 


r  H  dv 
Vsamp,  rf 

r  .  H  dv 

•'cav.  rf 


was  evaluated  using  the  Bessel  integral  expression 


/  r  tJ„  (X  d  r  -  ^  {tJ„  (X  c)f  +  <X  c)f}  -  J„  (X  c)  J^,  (X  c) 


and  the  Bessel  expansion  to  2nd  order  for  small  a. 
The  result,  accurate  for 

^^01  ^  ^  .  . 


p  +  a  +  R  +  s 


where 


'X  /  2\  /2Xrt,  Iv  aV 

o  I  /  a  I  A  01  / 

X  /  Is  /  I  D  / 

g/  \  /  \  o  / 
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The  equation  used  to  estimate  the  loading  due  to  Ge  is 

cal*  E  ^  dv 
ft  B  o  «Jcav.  rf 

'^l  (•  <yft  r  2  . 

I  Ge  E  f  dv 

t<samp.  rf 

where  °Ge  is  the  Ge  conductivity  and  is  assumed  independent  of  E  so  that  we  may 
remove  it  from  the  integral..  Further,  °^Ge  is  assumed  small  enough  so  as  not  to 

I 

affect  the  field  configuration.  The  calculation  is  similar  to  that  of  n  and  both 
results  are  plotted  in  Fig.  7.  For  =*  3.19  cm  (f  “  9.4  kmc)  we  have  D  a:  1 .20  cm. 
kg.  and  kg^  are  taken  to  be  11.8  and  16  respectively. 

The  unloaded  Q,  of  the  cavity  is  0^  —  8700  and  and  Q,j  are  added  in  parallel  to 
obtain  the  loaded  Q.  In  view  of  the  low  conductivity  of  intrinsic  Ge  at  room  temp, 
it  will  be  necessary  to  use  a  small  sample  and  low  temperatures  to  maintain  a  reason¬ 
able  0.  The  reduction  in  n,  of  course,  puts  a  lower  limit  on  sample  radius. 


I.  Nisenoff  and  Fan,  Semiconductor  Research,  Fifteenth  Quarterly  Report,  Feb.  1, 

1960  to  April  30,  1960. 

OPTICAL  BEHAVIOUR  OF  IRRADIATED  SILICON 
A.  K.  Ramdas  and  H.  Y.  Fan 

Our  previous  work^  has  shown  that  the  two  infrared  absorption  bands  at  11.98 
microns  and  11.56  microns,  produced  in  silicon  by  neutron  or  electron  irradiation, 
are  associated  with  the  presence  of  oxygen  in  the  material.  The  previous  study  of 
the  growth  of  these  absorption  bands  with  4.5  Mev  electron  irradiation  indicated  that 
the  11.98  micron  band  grew  at  about  the  same  rate  but  saturated  at  a  lower  level  with 
irradiation  at  liquid  nitrogen  temperature  as  compared  with  irradiation  at  0°K.  Such 
behavior  was  difficult  to  explain  on  a  simple  basis.  However,  the  low  temperature 
cell  used  in  the  experiment  was  not  suitable  for  accurate  optical  measurement,  and 
the  indication  for  the  saturation  behavior  was  given  by  only  one  point  in  the  curve. 

I 

The  Irradiation  at  liquid  nitrogen  temperature  has  been  repeated  with  a  low 
tenf>erature  cell  of  better  design.  As  was  found  In  the  previous  measurement,  the 

II. 98  micron  band  grew  at  first  linearly  with  electron  flux,  at  about  the  same  rate 
as  with  O^K  irradiation.  Furthermore,  the  band  did  not  saturate  at  a  lower  level  as 


1 
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indicated  in  the  previous  measurement  but  tended  to  saturate  in  the  same  manner  as 
with  0°K  saturation.  Thus  there  is  no  appreciable  difference  between  the  two 
temperatures  with  regard  to  the  production  of  the  11.98  micron  absorption  band  (See 
Fig.  8). 

in  agreement  with  the  previous  measurement,  the  11.56  micron  band  was  not 
observed  during  the  liquid  nitrogen  irradiation.  The  band  was  observed  after  the 
sample  had  been  warmed  up  to  room  temperature.  Irradiation  at  liqiuid  helium  temper¬ 
ature  is  planned  in  order  to  ascertain  whether  the  formation  of  the  11.98  micron  band 
is  inhibited  at  sufficiently  low  temperature. 


1.  A.  K..  Ramdas  and  H.  Y.  Fan  -  Semiconductor  Research,  First  (Quarterly  Report, 
October  1,  1960  to  December  31,  1960. 
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III.  LOW  TEMPERATURE  STUDIES 

SPECIFIC  HEAT  MEASUREMENTS  OF  INDIUM  WITH  SMALL  PERCENTAGES  OF  IMPURITY 
P.  H.  Keesom  and  B.  J„  L.  van  der  Hoeven,  Jr. 

Last  year  Bryant  and  one  of  us^  found  that  the  specific  heat  of  superconductive 
indium  between  0.4^K  and  0.7^K  was  less  than  the  lattice  contribution  In  the  normal 
state . 

2 

Ferrel  proposed  an  explanation  of  this  specific  heat  anomaly  which  is  based  on 

an  anomalous  dispersion  In  the  phonon  spectrum.  The  interaction  between  phonons  and 

electrons  leads  to  a  self  energy  of  the  phonons,  due  to  creation  and  reabsorption  of 

electron-hole  pairs.  The  energy  gap  In  the  electronic  states  of  the  superconductor 

prevents  the  creation  of  electron-hole  pairs  for  phonons  with  too  small  an  energy, 

and,  consequently,  the  self  energy  term  disappears  for  a  certain  range  of  phonon 

frequencies.  This  leads  to  a  lowering  of  the  frequencies  of  the  phonons  in  the 

superconductive  state,  relative  to  the  phonons  in  the  normal  state.  As  a  consequence, 

2 

a  new  term,  equal  to  -BT  ,  is  obtained  In  the  specific  heat  of  the  superconductor. 

A  distinct  advantage  of  this  theory  Is  that  It  predicts  the  effect  for  indium, 
but  not  for  tin,  since  it  depends  on  the  coupling  between  phonons  and  electrons. 

As  one  of  the  experimental  tests,  Ferrell  indicates  that  for  phonons  with  wave¬ 
lengths  longer  than  the  coherence  length,  this  effect  should  disappear.  The  coher¬ 
ence  length  Is  a  measure  of  the  extension  of  the  electron  wave  packet  taking  part  in 
the  superconductive  phenomena.  His  idea  is  that  for  phonon  wavelengths  longer  than 
the  coherence  length,  we  expect  the  normal  bulk  elastic  constants  to  give  a  correct 
description  of  the  sound  waves,  and  thus,  also  the  normal  lattice  specific  heat. 

3 

From  the  work  of  PIppard  ,  Ferrell  deduces  that  if  a  metal  (such  as  indium)  showing 
the  heat  anomaly  were  heavily  enough  doped  with  an  alloying  impurity  to  reduce 
significantly  the  coherence  length,  then  the  anomaly  would  disappear. 

Measurements  on  indium  doped  with  tin  are  now  in  progress. 


1.  C.  A.  Bryant  and  P.  H.  Keesom,  Phys.  Rev.  123.  No.  2,  491  (1961). 

2.  R.  A.  Ferrell,  Phys.  Rev.  letters  6^  No.  10,  541  (1961). 

3.  A.  B.  Pippard,  Phil.  Mag.  1104  (1955). 
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P lEZO-THERmL  CONDUCT IV ITY  IN  GERMANIUM 
R.  J.  Sladek 

1  y 

The  thermal  conductivity  of  degenerate  n-type  germanium  (containing  2  x  10  Sb 

atoms/cc)  has  been  measured  between  1 .5°K  and  4.1°K  with  tensile  stresses  from  zero 
8  2 

to  1.5  X.  10  dynes/cm  applied  in  the  [HO]  crystallographic  direction.  Fig.  9  shows 

the  temperature  dependence  of  the  thermal  conductivity  with  no  applied  stress  and 
8  2 

with  1 .0  X  10  dynes/cm  applied.  For  all  tenperatures,  this  applied  stress  increases 
the  thermai  conductivity  only  slightly  (about  lO^O.  This  is  in  sharp  contrast  to  the 
behavior  of  germanium  more  lightly  doped  with  Sb  for  which. large  increases  in  the 
thermal  conductivity  can  be  obtained  at  the  lower  temperatures'.  A  detailed  compar'- 
ison  of  these  effects  Is  given  in  Fig.  10  which  shows  that  the  ratio  of  unstressed 
to  stressed  thermal  conductivity  is  independent  of  temperature  for  the  degenerate 
sample  while  this  ratio  decreases  strongly  with  decreasing  temperature  for  the  more 
lightly  doped  sample  (data  on  the  latter  from  reference  1). 

The  relatively  small  piezo-thermal  conductivity  effect  in  the  degenerate  sample 
can  be  explained  qualitatively  by  noting  that  in  such  a  sample  the  donor  levels  have 
merged  with  the  conduction  band  and  the  donors  are  ionized  even  at  the  lowest 
temperatures.  Consequently,  the  strongly  stress -dependent  scattering  of  phonons  by 
bound  donor  electrons'  is  absent. 

The  different  temperature  dependences  of  the  piezo-thermal  conductivity  for  the 
degenerate  and  non-degenerate  samples  also  suggest  that  the  mechanisms  by  which 
electrons  scatter  phonons  are  different  in  these  two  types  of  materials. 


1.  R.  U.  Keyes  and  R.  J.  Sladek,  Phys.  Rev.  (Jan.  15,  1962). 


GERMANIUM  RESISTANCE  THERMOMETERS 
N.  Pearlman  and  H.  P.  Mathur 

In  the  measurements  described  in  earlier  Progress  Reports,  it  was  found  that  the 
reproducibility  of  the  thermometers  in  successive  runs  was  rather  good,  provided  that 
relative  values  R  (T)/R^  were  compared,  where  R^  is  the  resistance  at  an  arbitrary 
temperature  (usually  R^  “  R  (4°K)).  However,  R  (T)  itself  was  usually  reproducible 
only  to  such  an  extent  that  a  given  resistance  value  might  correspond,  on  different 
runs,  to  temperatures  differing  by  as  much  as  A  T  ~  0.l**K. 
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Part  of  these  differences  could  be  accounted  for,  on  runs  between  which  leads 
had  been  resoldered  to  the  thermometer,  by  the  possibility  that  melting  of  the  solder 
attaching  the  heavy  leads  to  the  thermometer  had  changed  the  effective  sample  dimen¬ 
sions  between  leads.  To  avoid  this,  the  heavy  leads  were  attached  to  the  thermometer 
with  tin-lead  solder  while  the  thin  wires  were  attached  to  the  heavy  leads  with  low- 
melting  cerroseal  solder. 

Another  possible  factor  affecting  reproducibil ity  could  be  strain  in  the  thermo^ 
meter  due  to  differential  thermal  contraction  between  the  thermometer  and  the  base  on 
which  It  is  mounted,  or  the  adhesive  bonding  It  to  the  base.  The  former  is  apparently 
not  too  important,  since  thermometers  mounted  on  a  6e  base  did  not  behave  differently 
from  those  measured  earlier,  mounted  on  a  copper  base,  in  order  to  eliminate  the 
latter  influence  thermometers  were  made  in  the  "L"  shape  shown  in  Fig.  II  and  mounted 
on  the  short  leg.  With  this  configuration,  reproducibi I ity  on  successive  runs  becomes 
a  few  mil  1 idegrees.  The  thermal  contact  resistance  dT/dP  between  thermometer  and 
base,  where  P  is  power  dissipated  in  the  thermometer,  is  about  the  same  as  that  for 
the  earlier  thermometers,  roughly  a  mlilidegree  per  microwatt,  which  is  about  two 
orders  of  magnitude  smaller  than  that  reported  for  encapsulated  Ge  thermometers  by 
Llndenfeld^ . 


1.  P.  Lindenfeld,  Rev.  Scl.  Inst.  32,  9  (1961). 
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IV.  GENERAL 


SHALLOW  ACCEPTOR  STATES  IN  GERMANIUM 
H.  M.  James  and  K.  S.  Mendel  son 


Solutions  of  the  effective  mass  equations  for  shallow  acceptors  in  germanium 
have  been  obtained  using  a  variational  method  in  which  the  angular  dependence  of  the 
trial  function  was  determined  from  symmetry  considerations'  and  the  radial  dependence 
was  determined  by  the  variational  calculation. 

Let  J  ,  J  ,  J  be  matrices  corresponding  to  components  of  angular  momentum  3/2. 

X  y  2 


fa  •’J  =  i  (ja  Jp  +  Jp  •>«)  ■ 


Then  the  effective  mass  Hamiltonian  for  acceptors  in  germanium  can  be  written 


where  Y, ,  Y^,  Yo  are  constants  which  can  be  determined  from  cyclotron  resonance  exper- 
2  3  ^  ^ 

iments  ’  .  The  energy  has  been  chosen  to  be  positive  in  the  direction  of  increasing 
hole  energy;  thus  the  bound  state  energies  of  the  acceptor  will  be  negative. 

The  ground  state  wave  functions  transform  according  to  the  Q  irreducible 
representation  of  the  double  tetrahedral  group.  The  function 


having  this  transformation  property,  was  used  as  a  trial  function'.  The  first  term 

has  s-llke  dependence  on  angle,  and  the  others  have  d-like  dependence.  (Factors  of 
2 

r  were  included  in  the  d-like  angular  terms  for  convenience  in  calculation.) 
Extremal izing  the  functional 
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L  -  ;  >),*  (H  -  E)  I  d\ 

with  respect  to  variations  of  the  f's  leads  to  a  system  of  ordinary  differential 
equations  for  the  f's.  It  is  convenient  to  define 


g,  -  rf,, 

_  3  , 

”  ''  ^2’ 
,  3  , 

93  r  f^. 


The  differential  equations  for  the  g's  are  then 


2/3 


gj"  +  (e  +  ^)  g,  -  5^^  (92"  +  ^  92*)  +  5^^ 


(g,"  +  7g3‘> 


2  n/S* 


^(g,''-7  3i'-^%5,)  +  g2"+(.  +  f-\)92 
M  r  r 


'1 


<91”  ■  7  9|’  %  9|’  93"  +  <•  +  7  ■  93 

r  r 


0, 

0, 


(3a) 

(3b) 

(3c) 


Here  we  have  used  as  units  of  energy  and  length 

4 


m  e 


2  2 
2  r  h  Y 


2  ’ 
K  h  Y 


m  e 


=  0.00403  ev. 


o 

no  A. 


The  functions  gj  roust  satisfy  the  boundary  conditions 

g,  -•  0  as  r  -»  0,  r  •  (I  *  1,  2,  3). 


We  first  observe  that  these  equations  have  solutions  with  g^ 
satisfying 


0  and  g^,  g^ 


92”  +  («  +  7  “  %)  9 


r  2'  ^2 
r 


h  “  “T^yT  ®2- 


0, 
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This  equation  for  92  Is  simply  the  equation  for  hydrogenic  radial  functions  with 
1=2.  It  will  have  acceptable  solutions  for  e  *  -1/n^,  n  *  3,  4,  5,  ...  These 
solutions,  however,  do  not  give  an  approxinuition  to  the  ground  state  function.  The 
lowest  energy,  «  “  -1/9,  E  •  -0.00045  ev,  is  much  higher  than  the  experimental  ground 
state  energies^  of  -0.010  to  -0.014  ev. 

If  ^  0,  we  can  eliminate  it  from  Eqs.  (3a)  and  (3b).  This  gives 
2 

^~2  (  Y3  92  ^2  93)  +  («  +  ■^  -  ■^)  (  T3  §2  ^2  93)  “  0. 

dr^  r 


This  is  again  the  hydrogenic  radial  equation.  There  will  be  acceptable,  non-trivial 
solutions  for  ^3  Y3  ^  ^2  ®3  when  c  "  -1/n^.  As  we  have  seen,  these  solutions 

will  not  be  approximations  to  the  ground  state  functions.  Thus  for  the  ground  state 
we  must  have 


93  “ 


92- 


Elimination  of  g^  from  Eqs.  (3)  gives 


g,"  +  (*  +  7)  9,  -  A  (92"  +  7  92') 


B  (g,"  -  7  9,'  +  %  9,)  +  92"  +  <•  +  7  -  \>  92  “  0. 


where 


2  (2y2^  +  3Y3^) 

5  Yi  Y2  ' 


(4a) 

(4b) 


B 


Yi 


These  equations  cannot  be  solved  analytically.  We  have  therefore  obtained  numerical 
solutions  using  the  Runge-Kutta  method^. 

To  determine  the  initial  conditions,  one  can  try  solutions  of  the  form 


9l 


92 


Si  a  r**^, 

T-0 


E  b  r^ 


(5) 
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One  finds  that  n^  =  n2  “  n  3,  1,  0,  -2.  We  note  that  the  indices  differ  by 
integers.  In  this  case,  there  will  not  be  a  solution  of  the  form  (5)  for  each  value 
of  n^.  There  will  be  a  solution  of  the  form  (5)  for  n  ^  3.  We  denote  this  solution 
by 


Yl  “ 

In  addition  there  will  be  a  solution  of 

H  “ 


<r*0 


E  b 


T“0 

the  form^ 
y,  In  r  + 
72  1”  r  + 


91 

92 


(6) 


We  denote  this  solution  by 


z 


I 


=^2 


y,  in  r  +  r  S  c  r**^, 

<T“0  ^ 

2  *  T 

y,  in  r  +  r  Ed  r^. 

^  T-0  ^ 


(7) 


From  the  differential  equations,  one  can  derive  recurrence  relations  for  the  coeffi¬ 
cients  a  ,  b  ,  c  ,  d  .  There  are  two  other  linearly  independent  solutions  of  Eq.  (4), 
T  T  T  T 

but  they  do  not  behave  acceptably  as  r  0. 

It  is  also  necessary  to  consider  the  asymptotic  behavior  as  r  ■  There  are 
four  linearly  independent  asymptotic  solutions  that  have  the  forms 


or 


+  ar 

g,  ~  aer-  , 


g-  ~  be- 


+  Ctr 


.be;fcp^ 


(a  >  P) 


(8) 


Op 

As  r  -*  ®,  the  solution  with  e  will  blow  up  more  rapidly  than  that  with  e  hence 
this  will  be  called  the  dominant  solution,  and  that  with  e*^  the  sub-*dominant 
solution. 

For  a  given  energy,  one  can  use  the  solutions  described  by  Eqs.  (6)  or  {7),  or 
any  linear  combinations  of  them,  to  determine  gj(r^),  9*|Cr^), 
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(r  small)  for  a  solution  that  Is  well-behaved  as  r  -*  0.  Starting  from  r  with  these 
o  o 

Initial  values,  one  can  Integrate  Eq.  (4)  numerical iy  to  obtain  a  solution  that  Is 
well-behaved  at  the  origin  but  will,  in  general,  become  infinite  as  r  -*  «».  If  an 
arbitrary  1  inear  combination  of  Eqs.  (6)  and  (7)  Is  used  to  obtain  the  initial  condi¬ 
tions,  the  solution  at  large  r  wi 1 1  approach  some  linear  combination  of  the  dominant 
and  sub-dominant  solutions.  In  general,  such  a  solution  will  blow  up  even  if  e  is  an 
eigenvalue  of  the  problem.  To  obtain  a  solution  that  is  well-behaved  as  r  -*  ■>,  it 
is  necessary  to  adjust  not  only  the  energy,  but  also  the  linear  combination  of  Eqs. 

(6)  and  (7)  that  is  used  to  obtain  the  Initial  conditions. 

One  can,  for  any  t,  choose  the  initial  conditions  so  that  only  the  sub-dominant 
term  appears  at  large  r.  For  an  energy  that  Is  an  eigenvalue,  the  sub-dominant  term 
disappears  also,  and  one  obtains  a  solution  that  is  well-behaved  as  r  To 

eliminate  the  dominant  solution,  the  following  procedure  was  used.  For  each  energy 
two  solutions  (g|,  g2^  were  determined,  starting  from  initial  conditions  determined 
by  Eqs.  (6)  and  (7),  respectively.  It  was  found  that  the  solution  given  by  Eq.  (6) 
had  essentially  pure  dominant  behavior  at  large  r.  We  took  the  ratio  9^/92 
solution  at  large  r  to  be  a/b  (Eq.  (8)).  Using  the  value  of  a/b  obtained  in  this  way, 
we  chose  a  linear  combination  of  the  solutions  such  that  the  ratio  9^/92 
of  the  range  of  integration  (determined  by  practical  considerations  of  computing  time) 
was  -a/b.  The  resulting  curves  gave  essentially  pure  sub-dominant  behavior  at  large 
r.  The  value  of  gj  at  a  point  near  the  end  of  the  range  of  integration,  chosen  so 
far  out  that  a  well-behaved  solution  would  be  negligible,  was  used  as  a  measure  of 
the  amount  of  sub-dominant  solution  present.  By  repeating  this  process  for  several 
energies,  it  was  possible  to  find  eigenvalues  and  corresponding  eigenfunctions. 

Using  this  method,  we  determined  the  energy  and  radial  functions  for  the  ground 
state.,  and  for  the  first  excited  state  described  by  the  trial  function  (2).  We  also 

I 

applied  a  similar  method  to  the  and  excited  states  .  For  the  states,  we 
used  the  trial  function 
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-  I  C7(z^  -  r^)(x  -  ty)  -  (x  +  ly)^] 

31z  (s^  -  I  rh 

I  c<z"  -  i  r^Xx  +  iy)  +  (X  -  ly)’] 

tz(x  +  ly)^  +  z(x  -  ly)^] 

Here  there  are  only  two  terms  In  each  trial  function,  and  the  variational  procedure 
gave  immead lately  two  coupled  differential  equations  for  the  f*s.  These  equations 
were  solved  numerically  by  the  method  described  above.  The  energies  thus  found  are 
shown  In  Table  I,  and  the  radial  functions  in  Figs.  12  -  16,  These  results  were 
obtained  using  the  values  of  »  ^2'  ^3  9*^®”  reference  3. 


3  (x  -  ly) 
2  Iz 
X  +  iy 
0 


f^  (r)  + 


TABLE. I .  Eigenvalues  of  Effective  Mass  Hamiltonian  for  S ha I  low  Acceptor  States  In 
Germanium. 


Trial  Function 


r, 

17 


Energy 


•o 

o 

1 

(D 

< 

''-2.32 

-9.28 

-0.690 

-2.76 

-0.484 

-1 .94 

,-0.250 

-1 .00 

-0.146 

-0.584 

1.  D.  Schechter,  Theory  of  Sha I  low  Acceptor  States  in  Silicon  and  Germanium.  Ph..D. 
Thesis,  Carnegie  Institute  of  Technology,  unpublished. 

2.  G.  Dresseihaus,  A,  F.  Kip,  and  C.  Kittel ,  Phys.  Rev.  3.68  (1955). 

3.  B.  W.  Levinger  and  D.  R.  Franki,  J.  Phys.  Chem.  Solids,  20,  281  (1961). 

4.  P.  Fisher  and  H.  Y.  Fan,  Phys.  Rev.  Letters  2,  456  (1959). 

5.  F.  B.  Hildebrand, "Introduction  to  Numerical  Analysis^  (McGraw  Hi  1 1 ,  New  York, 

1956)  pp.  233-239. 

6.  E.  L.  Ince,  "Ordinary  Differential  Equations"  (Dover  Publications  Inc.,  New  York., 
1956)  p.  160  and  Sec.  15.311. 
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GROWTH  OF  Z  INC  TELLUR  IDE  SINGLE  CRYSTALS  FROM  THE  VAPOUR  PHASE 

P.  H.  Klose 

Work  has  been  started  to  prepare  the  compound  semiconductor  zinc  teiiuride,  ZnTe.. 

Equimolecular  parts  of  the  chemically  pure  constituents  are  weighed,  mechanicaily 

mixed,  and  meited  together  in  a  quartz  tube.  A  purified  hydrogen  ambient  is  used  to 

prevent  oxidation  during  this  process.  Any  excess  metal  is  distil ied  away  afterwards 

in  a  high  vacuum  furnace.  The  material  so  prepared  is  already  quite  dense  and  is  In 

the  shape  of  a  rod.  This  rod  is  then  placed  in  a  furnace  for  the  growing  of  the 

crystal.  The  furnace  is  essentially  similar  to  that  described  by  Polich  and  co- 

workers^.  Experiments  are  being  carried  out  to  determine  the  optimum  temperature 

profile  of  the  growing  furnace  and  the  corresponding  rate  of  growth  in  order  to  obtain 

single  crystals.  Measurements  of  Hail  coefficient  and  resistivity  of  the  material 

2 

prepared  have  been  made  at  room  temperature.  Hall  mobilities  of  50  cm  /volt-sec  have 
been  found  for  samples  with  resistivities  of  the  order  of  I  0  cm. 

i.  Piper,  Polich,  J.  A.  P.,  3^,  1278  (1961). 

MATERIAL  PREPARATION 
K.  Masumoto  and  L.  Roth 

Work  has  continued  on  the  preparation  of  large  crystals  of  ZnSb  under  various 
growth  conditions  using  the  Bridgman  method.  Attempts  were  made  to  obtain  single 
crystals  by  the  Czochraiski  technique  using  the  Bridgman  method  crystals  as  seeds. 

So  far,  it  has  not  been  possible  to  obtain  large  single  crystals  without  twinning. 
Improvements  are  being  made  in  the  crystal  pulling  apparatus  with  the  aim  of  obtaining 
single  crystals. 

The  usual  work  of  making  single  crystals  of  Ge,  InSb,  and  GaSb  continues. 
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